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O ABSTRACT 

When the gravitational lensing potential can be approximated by that of a circularly symmetric 
system affected by weak perturbations, it is found that the shape of the resulting (tangential) 
£h , caustics is entirely specified by the local azimuthal behaviour of the affecting perturbations. 

This provides a common mathematical groundwork for understanding problems such as the 
close-wide (d «-» d~ l ) separation degeneracy of binary lens microlensing lightcurves and the 
shear-ellipticity degeneracy of quadruple image lens modelling. 
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O 1 INTRODUCTION 
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The recent announcement bv lBond et alj 12004) of the most convincing planetary micro lensing event iMao & Paczvnskill99ll) to date is one 
of the greatest success stories of the nearly two-decade old promise of microlensing IPaczvnsklll98d) . However, despite such successes, 
the study of microlensing i s still hampered by various hurdles, one of the most persistent being the problem of degenerate lightcurves (e.g., 
lDominilJl999J : lAfonso et al 1200(1) - that is, several distinct physical systems can be used to model observed lightcurves. 

For the purpose of formal discussion, the degeneracy of microlensing lightcurves may be categorized into two separate but related 



i 

phenomena. The first, which may be termed as 'external' degeneracy, is caused by imperfect observations yielding highly correlated mea- 
surements of the parameters that control the observed lightcurves. For instance, in most highly blended events, without extremely good 
coverage of certain wing portions of the lightcurve, one only expects to measure the ratio of the time-scale to the peak magnification (or 
the product of the time-scale and the source flux) but not the parameters individually (c.f.. lGoul(il99tih . On the other hand, the second type 
of degeneracy, which will henceforth be referred to as 'internal' degeneracy, is due to the fact that the number of independent lightcurve- 
controlling parameters is smaller than the number of parameters that are used to describe the underlying lens systems. In other words, 
j-^ . even with perfectly specified lightcurves (consequently, perfectly determined lightcurve parameters), the underlying lens system may not 
be uniquely recovered. This internal degeneracy of microlensing lightcurves can be further divided into two broad classes; the 'intrinsic' 
degeneracy among magnification structures from different systems, and different lens systems tracing identical or intrinsically degenerate 
paths in the given magnification structure (hereafter the 'extrinsic' degeneracy). The most widely acknowledged degeneracy of microlensing 
lightcurves, that is, the 'time-scale degeneracy' relating to the mass of the lens, the distances to the source and the lens, and the relative 
motions between the observer, the source, and the lens , is an exam ple of the e xtrinsic dege neracy. Other examples of the extrinsic degeneracy 
are the discrete parallax degeneracies fe.g.. lSmith. Mao."& Paczvnskill2003t iGouldl 2004). However, it is the intrinsic degeneracy that is of 
the most theoretical interest as it can be studied through the analysis of the lens equation as a whole without specified reference to the path 
of the system through the magnification structure or the observational conditions. Furthermore, understanding of the intrinsic degeneracy 
can greatly facilitate the identification of various extrinsic and external degeneracies. The most obvious example of an intrinsic degeneracy 
is the azimuthal symmetry as well as the radial similarity of the lensing magnification of the point mass lens. The well-known result of 
the 'inverse-square-root' behaviour of magnification 'inside' fold caustics ("e.g.. lGaudi & Pettersl2002l) is a less obvious example of a local 
intrinsic degeneracy. 

It is notable that the intrinsic degeneracy is intimately related to the symmetry of the lens equation. In general, if (the potential of) the 
system possesses a certain symmetry, one can expect that there should exist an intrinsic degeneracy of magnification related to it. Note that 
the magnification is basically the second derivative of the potential. This also means that the study of the intrinsic degeneracy beyond the 
incidental case studies can be greatly systematized by recognizing the associated symmetry structure of the lens system. However, beyond 
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the obvious example of the circularly symmetric lens, studies of degeneracies rooted in symmetries have been minimal, at best. This paper 
is one of the first attempts to understand certain intrinsic degeneracies of microlensing by searching for a symmetry in the lens system. It 
is found that this approach proves itself by providing a unified scheme of understanding the problem of a certain well-known microlensing 
degeneracy and a seemingly remote problem of lens modelling of a quadruple images. 

The primary focus of this manuscript is a point mass lens (and a more general circularly symmetric lens in later sections) that is under 
the external influence of certain perturbations. It should be noted not only that the symmetry in such systems can be analysed with relative 
ease, but also that they can be applied to various realistic systems. In Section|2| the basic properties of the point mass lens are recapitulated, 
and in Section[3] a perturbative approach to gravitational lensing is developed. In particular, in section l3~Tl the perturbative approach is used 
to derive an approximate expression for caustics and critical curves. Following this, in Section|4] it is shown that there exists a certain set 
of perturbations that yields invariant caustics, which may be seen as the main finding of this monograph. This finding is applied to interpret 
one of the well-known microlensing degeneracy problems in the next Section [5] The restrictions imposed on the studied lens system are 
subsequently relaxed in Section|6| For example, the discussion is extended from a point mass lens under perturbations to a general circularly 
symmetric lens system under perturbations in Section l6^2l and from external perturbations to perturbations associated with the system itself 
in Section lo31 These two generalizations lay the basis for a study of the shear and ellip ticity in lens mo delling late r in the same section. 

In this paper, the lens equation formulated in terms of complex numbers iBourassa & Kantowski|[l975lwitJl990l : IWitt & Mao|ll995l) 
is extensively used. In Appendix IaI most of the mathematical terminologies regarding complex analysis found in this manuscript are sum- 
marized, and in Appendix IbI the basic gravitational lens theory is redeveloped using complex-number notation. The treatments given there 
are minimal and deliberately casual. More in-depth treatments of complex analysis includin g rigorous d efinitions and proofs can be found in 
any standard text of complex analysis (e.g.. lAhlforsll979T) or of mathematical methods (e.g.. lArfken & Webej200(il) . Interested readers may 
also consult references regarding potential theory in a plane. 



2 SCHWARZSCHILD LENS 

When the lens system can be approximated by a Schwarzschild metric, one can write the lens equation in complex notation (see AppendixlBl 

as 

4=Z-r, (2.1) 

z 

which relates the apparent angular position of the lensed image z to the (would-be) angular position of the source f s in the absence of the 
lens. Note that, throughout this paper, the complex conjugate is represented with an overline (or upper bar) so that z is the complex conjugate 
of z. Here, the lens is located at the coordinate origin, and all the angular measurements are made in units of the angular Einstein ring radius, 



<2R S 



i/2 



where R Sl . h = 2Gmc~ 2 is Schwarzschild radius of the lens mass m, and D rA = D^D L D & is the relative parallax distance 1 , and D s and D L are 
the distances to the source and to the lens respectively, while D LS is the distance to the source from the lens. 

It is relatively straightforward, if not entirely trivial, to solve the lens equation <2. 1> to find the image positions z for a given source 
position f s . However, one can discover a few interesting properties of the lens equation 12. H before actually solving it. For example, from 
the symmetry of the equation, one can immediately find that, if z\ is an image position for a given source position, zi = -z\ l should also be 
an image position for the same source position (also note that z^ = -Zi). In fact, it is easy to show that equation <2.H always allows only 
two images for any given source position ( s unless £, = 0, in which case the equation 12. H is reduced to the equation of a unit circle, |z| = 1. 
That is, the image becomes a ring (which is also naturally expected from the intrinsic circular symmetry of the system) with its radius given 
by equation I2.2t . This ringed image is sometimes referred to as 'Einstein ring.' 

One important quantitative measure in gravitational lensing is the magnification factor of lensed images. In a purely mathematical sense, 
the lens equation defines a mapping from the image position to the source position, both of which are vectors in two dimensional spaces. 
The local - differential - behaviour of this mapping therefore can be studied from its Jacobian, which is basically a linear transformation 
that approximates the lens mapping locally. It follows naturally that the Jacobian determinant of any given lens mapping is the inverse 
magnification factor of the lensed image in the limit of a point source. If the lens mapping is given by the lens equation <2.1> . one can easily 
show that the associated Jacobian determinant is that 

J s = Ot.or-. - (U\(l :-. = |3^ s | 2 - |9^ t | 2 = 1 - iy. (2.3) 

kl 4 

Here, 9 ; = (d/dz) is simplified notation for the partial derivative operator. Equation 12.31 gives the inverse magnification for a single image. 
If one wants to find the total magnification, accounting for all images for a given source position, one needs to find the sum of the inverse 
absolute values of the Jacobian determinants corresponding to all the images. It is also notable that the Jacobian determinant given by 



IfDLS = D s 



d- s K 
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equation 12.31 vanishes when \z\ = 1 (consequently f s = 0), which implies that the magnification for this image is formally infinite - the 
lensed image of a '0' -dimensional point source is a T -dimensional ring. In general, the loci in the image space defined by = are referred 
to as critical curves, and the image of the critical curves under the lens mapping in the source plane as caustics. Thus, for the point mass lens 
case, the caustic is a point (£, = 0), and the critical curve is a circle (|z| = 1). 



3 SLOWLY- VARYING SMALL NULL-CONVERGENT PERTURBATION 

Let us consider the situation in which a point mass lens is subject to some null-convergent - i.e., the continuous surface lens mass density 
is zero except finite number of isolated mass points - perturbations. Then, the perturbation part of the lensing potential if/ satisfies the two- 
dimensional Laplace equation V 2 tf/ = 0. Since the real and the imaginary parts of any complex analytic function are harmonic, there exists 
a complex analytic function i// c (z) whose real part is the same as tf/ and whose imaginary part is a solution of first-order partial differential 
equations derived from the Cauchy-Riemann condition (see Appendix lAl. Since 2if/(x, y) = tf/ c (z) + i[>c(z) and 29=t/' = d x if/ + id y i/r (see 
Appendix[B) where z = x + iy, the resulting lens equation can be written in complex notation as 

S = £,-[ej& + c], (3.1) 

using a complex analytic function e/(z) = ifr' c (z)—c, where c is a complex constant (without the loss of generality, e may be taken as real). Note 
that here and throughout this paper, the use of primed symbols is exclusively reserved either for the complex total derivative of an analytic 
function or for the ordinary derivative of a real-valued single-real- variable function, with respect to their argument, while the argument will 
be dropped whenever there is little ambiguity. 

If |e/(z)| <k 1 for |z - zol <k 1, where zo is the image position of a point mass lens satisfying ( + c = zo - Zg 1 , one can find the image 
displacement Sz = z - Zo caused by the small perturbation of ef(z) for a given source position, by a series expansion of equation 13. 1 1 up to 
first order of both Sz and e, 

8( = a^lofe + d&loSz - e/tfo) = 8z+%- efo, (3.2) 



where f = f(zo) = f(za). If J7j,o = 1 - kol 4 ^ (i.e., kol + 1), equation J3.2I is invertible for Sz. Setting S( = 0, this inversion leads to 

&=7F-U-§1 (3-3) 
Js.O \ z l 

That is, the new image forms where the perturbation is counterbalanced by the image of the local linear mapping that approximates the 
(unperturbed) lens mapping. 

The Jacobian determinant of the lens mapping J3.1> is 

J=l ~w +e ii + iY £2ff =Js+ w ^ f + ff ) ~ e2|/f ' (3 - 4) 

where J~ s is the part of Jacobian determinant that maintains the same form I2.3i as the point mass lens. The change of the value of the 
Jacobian determinant with the perturbation is therefore understood as the sum of two contributions (up to the linear order); one due to the 
direct additional contribution from the perturbation (e-term in eq. I3.4> and the other due to the change of value of J~ s at the new image 
position 

S(J S ) = d : J s \ Sz + d- z J s \ Sz =All- + ^]- (3-5) 



kol 4 \Z0 Zo, 

When the small-perturbation solution (eq. l3.3> is valid, one may further substitute the solution for Sz', 

c ( rr,_ 2e (fa i fo\_ le Zofo+Zofo 

Us> (l + koP)kol 2 Uo zo) kol 4 1 + kol 2 • { ' 

For this case, the final first order change of the Jacobian determinant with respect to the point mass lens case is given by the sum of the two 
terms 

2e / Zq/q + zlf ' zo/o + Zp/ o ^ 
kol 4 I 2 1 + ko| 2 

where zo is the position of the unperturbed image. In addition, if \6$\ <k LTs.ol = |1 - kol 4 L then the first order change of the magnification 
can be found by 

8 J = 2pe I zg/ ' + zgio' zo/o + zo/ o ) 

J Jl 2 - (kol 4 + koh 4 ) I 2 1 + kol 2 

where p is the parity of the image. Provided that the perturbation does not change the parity of the image, p = 1 if kol > 1 and p = -1 if 
kol < 1. Note that the effect from equation 13.81 is most likely to be negligible in practice. This is because 6 A ~ e, and furthermore, for both 
of the most relevant regimes of equation 13.81 . ^'({"'M) — > as zo — > 0, and zo(^ [ SA) — > as Zq 1 — > 0. 



8J = J-Js,o=^\^^ + ^-^), (3.7) 



r. T ' MV0 , <-VJU 1 «J« f Q ON 

5A = -ptpt = o 7T. ,4 , Z i-^ I ^ + i , ,_ ,9 1. (3.8) 
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3.1 caustic and critical curve 

As |zo| — > 1 and consequently J s0 — > 0, the perturbative solution (eq. l3.3t grows and eventually diverges. That is, the linear perturbation 
becomes invalid or incomplete when the image for the point mass lens approaches the Einstein ring (= the critical curve). This is because the 
local linear mapping (eq. l3.2> that approximates the lens mapping of the point mass lens becomes projective at the critical point. This implies 
that the general perturbative solutions at an arbitrary critical point require consideration of higher order effects. However, for the limited case 
in which the direction of the source displacement coincides with the projective axis, the solution can be obtained from the linear effect alone. 
Moreover, for the point mass lens, the caustic point is multiply degenerate, and therefore, for source positions near the lens position (= the 
caustic point), one can select a valid base point for the series expansion among any of the critical points on the Einstein ring that result in the 
source displacement along the projective axis. 

Suppose that the lens equation 13 . 1 i is series-expanded at zo = e 1 *- Then, with Sz = (Sr + iS(/>)e" 1 ', the linearized lens equation becomes 

<T = (Sr + i(ty)e'* + i ef(e-"") = 2<5re'* - efieT^), (3.9) 

where f = S£ because zo - Zq 1 =0. Since Sr e R, the linearized lens equation 13.91 has solutions if e~"*[£ + e/(e~"*)] 6 R or equivalently t/> is 
the solution of 

m = C~ + e [fa - *) - e 2 ' ^)] =0. (3.10) 

In general, equation 13. 101 allows multiple solutions in [0, 2n). In principle, for any given expansion base point zo = e"\ one can recover 
multiple solutions of Sz corresponding to each of the solutions of equation 13. 1UI if higher order effects are considered. 
From equations I3.4i and 13.51 . the Jacobian determinant corresponding to equation I3.9i is 



- + * + | (zo/o + &o) = 4<5r + e [e 2 '*/'(e^) + e^/V*)] = 2e* [( - f e («] , (3.11) 

Zo Zo 2 v ' 1 1 



J kol 4 
where 

CM) = ~\ [e 3l V(e i0 ) + e^/V*) + 2^)] . (3.12) 

Here, J7j,o = 1 ~ kol 4 = so that J = SJ, and also e R if <p is the solution of equation 13. 10i . Note that $ < 1 for any null-convergent 
lens system so that equation 13. Ill is only valid if |£ - ( c (<f>)\ <k 1 /2. The total magnification for the given source position can be obtained by 
adding the inverse of all Jacobian determinants corresponding to each solution of equation 13. 101 in [0, 27t), 

MO= J] " • (3 ' 13) 

*e(0^<27T|fW=0} ~ CcWJI 

Now note that if there exists t/> e [0, 27t) such that £ c (<f>) = C, then S(0 = and therefore A(0 diverges. In other words, (^(0)10 6 [0, 27t)j 
defines (the linear approximation of) the caustics - that is, £•(</>) is the parametric form of the (linear approximation of the) caustics. The 
(parameter for) cusp points formed along the caustics can be found by solving d(Jd<f> = for the parameter <f>. Since 

^ = - | [3ie 3i <7'(e i < 4 ) + ie^/V*) - ie^/V*) - ie- 2i <*/"(e^) - He-'* fie-'" 1 ')] 
= £El {[3e 2i V(e 10 ) + e^/V*)] - [3e _2i */'(e- ¥ ) + e _3i *f (e""'*)]} 

= ee i(4 3[3e 2i(4 /'(e 1(4 ) + e 3i */V*)] (3. 14) 

where 3[z] '■ C — > R is the imaginary part of z, the condition for the parameter (f> to define a cusp point is that 3e 2l */'(e'*) + e 3l(4 /"(e 1(4 ) 6 R. 

The expression for the corresponding (linear approximation of the) critical curve can be easily found by setting J = in equation 13. 1 1 i . 
that is, 

Sr = e X \ (0) = -£ [f&fXt!*) + e^/V*)] . (3.15) 

Strictly speaking, this only defines the local linear approximation of the critical curve (near zo = e 1(4 )- However, one may simply consider 
( 1 + <$r)e'* as a global expression for the critical curve parametrized by 0, 

zM = e*[l + e X , (</>)] = ^-\ [e 3,(S /'(e'*) + e^/V*)] • (3-16) 
Then, the image of z c (0) under the linearized lens mapping 

^.(0) = 2e'^ Xl (0)-f./le- i *), (3.17) 
indeed recovers the previous derived expression 13. 1 21 for the caustic. 
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4 INVERSE SYMMETRY OF PERTURBATION POTENTIAL 
4.1 caustic invariant perturbation pair 

If f(z) can be represented by a complex monomial az" near the unit circle \z\ = 1 and \ea\ <k 1, its caustic can be found from equation J3.12I 
by 



- -U<P) = «ae 3 "V"- lw + nae^e (1 -" w + 2a& mip = rcae ( " +2W + (n + 2)ae 



(4.1) 



The resulting expression is a linear combination of two complex exponentials. Motivated by the relationship between the two exponents and 
the scalar coefficients of each exponential, let us consider the substitution m = -(n + 2). Then, 



p ("+2)i(t 



+ (n + 2)ae 



(m + 2)(-a)e-""* + m(-a)e' 



(m+2)i0 



(4.2) 



that is, equation <4. H is invariant under the exchange of n «-» —(n + 2) accompanied by a «-» -a. In other words, up to the linear approximation, 



the shapes of the caustics due to the perturbations f(z) = az" and f(z) 



-az 



-(n+2) 



are identical. Furthermore, the expression J3. 121 is 



linear to terms involving /, and therefore, the superposition principle implies that the caustic resulting from the polynomial perturbation 
ft(z) = Y,n a nZ n is the same as the one from another polynomial perturbation f p (z) = Y, n (—a„)z ~ ( " +2) . Moreover, one can easily show that 
equation 13.101 is also the same for both cases, and therefore equation <3.13> implies that this caustic invariance actually extends to the 
correspondence of the magnification associated with any source position in the vicinity of the caustics. 

In fact, the pair of 'linear caustic invariant perturbations' (hereafter, LCIP) is also connected by the symmetry between the corresponding 
lensing potential. To see this, let us consider two complex potentials with inverse symmetry; if/' c (z) and i//%(z) = ^(z -1 ). If ift' c (z) has a 
convergent Laurent series expression in a domain containing the neighbourhood of the unit circle - since tft' c (z) is an analytic function of z for 
any null-convergent perturbation, this condition is basically the restriction to the absence of poles in the neighbourhood of the unit circle -, 
then the perturbations due to these potentials are 



#fe) = ^ b mZ m -» f,(z) = #'fc> = mb ^' X = X a " z "' 



(4.3) 



(4.4) 



so that they are indeed an LCIP pair. In terms of the real potentials, the symmetry between the LCIP pair is more straightforward. First, 
a general (real) solu tion of the two-dimen sional Laplace equation - allowing singularities at either the origin or infinity - in plane polar 
coordinates (e.g.. lCourant & HilberJl962T) is, by harmonic expansion, 



V 



ao 
2 

-lnzz + b + -)_ j 



In r + Z? + ^ {^n r " H — 7" j cos n $ + {bn 1 * 1 H — ~ j s i n n< P 

b»> + *, + ± f] \( a ,y + («r* + **) + i (ty + H ( e 

n= 1 L \ / 



(a„ - \b n ye** + — 



a_„ + ib-„ 



+ (a n + ib„)r"e 



a- n - ib-„ 



where all the constant coefficients, a„'s and b„'s are real. The corresponding complex potential is easily found as 

\jf c = Oolnz + 2b + X (c n z" + c_„z~") , 

and therefore, the complex potential of its LCIP pair is 

V c =a lnz + 2b + Y i (c„ Z -" + c_„z") , 

where c„ = a„ + ib„. Finally, the corresponding real potential of the LCIP pair is 

V 



(4.5) 



(4.6) 



(4.7) 



-\nzz + b Q + -^ 



+ (a-„ - ife_„)r"e""* + — — — + (a_„ + ib-„)r"e 



| In r 2 + b, + 1 £ [(£ + a_„r") (e^ + e-*) + i + i_„r») (e^ - e*) 

n=l L \ / 

a In r + Z? + ^ + a-„r"j cos n<f) + + b-„r" 



(4.8) 

In other words, the perturbing parts of the real potentials of the LCIP pair are related to each other by an inverse symmetry with respect to 
the unit circle (i.e., Einstein ring); that is, 6if/ p (r,<f>) = (5i/r'(r~', <f>) in plane polar coordinates (r, if)), or 6t// p (r) = Sift'ir/r 2 ) in vector notation. 
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Harmonic expansion of the potential, in addition, demonstrates that the presence of an LCIP pair is related to the intrinsic symmetry of two 
dimensional harmonic function, that is, the azimuthal structure is invariant under the radial distance inversion. 



4.2 caustic invariance and magnification correspondence 

Following the discussion in the previous section, the natural question arises whether the magnification correspondence between the pair of 
potentials with the inverse symmetry extends to sources lying far from the caustic. For the zeroth order, the answer is affirmative because the 
magnification of sources lying far from the caustic is just small a perturbation on the point mass lens case. However, to examine whether the 
correspondence actually extends to the first order perturbations, one needs to consider the approximations of equations 13.21 . 13.71 . and 13.81 
that provide one with the lowest order non-trivial effects due to perturbations for these cases. 

From equation J3.8I . the magnification change associated with the perturbation f,(z) = az" is found by 



On the other hand, the magnification associated with the corresponding perturbation f p (z) 

pe I + 2 \ / _a 

2 - (Izol 4 + Izol- 4 ) T l + kor 2 jU + 



-r, 7 -("+2) - 



SA„ = '- \n+ — T + — - . (4.10) 

P "i /l,,|4 , |,„-4\ 1 . u„l-2 }\=n+i I v ' 



Next, one notes that, for a given source position, the symmetry of the point mass lens implies that the positive and negative parity images 
should be related by zo <-> -Zg 1 ■ Hence, the linear changes in magnification on alternative parity images by a pair of perturbations f,(z) 
and f p {z) are symmetric for even n - an odd potential - while they are antisymmetric for odd n - an even potential. In other words, the 
pair of perturbations f,(z) = az" and f p (z) = (-l)" +1 az~ ( ' 1+2) yields correspondent linearly perturbed magnifications for the source position 
far from the caustics. Using the same argument in the previous section, one may generally conclude that this pair is in fact related, in 
terms of potentials, with a property that i/f p (r) = -if/'i-r/r 2 ) or i//%(z) = — z~ ). In addition, one may also state that the magnification 
correspondence between the LCIP pair of odd potentials is stronger than that of even potentials. 



5 EXAMPLES 

5.1 planetary perturbations 

IBozzal Il999h showed the perturbative analysis provides an ideal method to study microlensing from planetary perturbations. Here, some 
well-known degeneracies of planetary microlensing are reexamined using the perturbative approach. In general, microlensing by a star with 
a planet can be described by the lens equation, 

f = z-4-^r. (5.D 

z z - Z p 

Here, the position of star is chosen as the coordinate origin, and Einstein ring radius corresponding to the stellar mass alone as the unit of 
angular measurements. The projected angular location of the planet is given by z p , while q is the mass ratio of the planet to the star. From 
the comparison to equation 13. 1 i . if < q <k 1 (and |z — z p | » q), equation 15. 1 i can be considered as a point mass lens under a small 
null-convergent perturbation that comes from 

/(z) = (z-z p r 1 ; f'(z) = -(z-z v Y\ (5.2) 

which is analytic everywhere except z = z p where it forms a pole. Then, if z p is sufficiently far from the circle |z| = 1, equation 13. 121 provides 
a parametric representation for the linear approximation of the (central) caustics; 

2 _ e 3i * e _i * 2 

q ~ (e^-Zp) 2 + (e-* - z p ) 2 " e-^-V ( ' 

where if> e [0, 2n) is a parameter. To examine the degeneracies of planetary microlensing related to these caustics, it is helpful to further 
manipulate the last two terms in the equation 15.31 algebraically, that is, 

e-* 2 z„ 2 e i0 2z p 'e'* ( T' Y .. e 3 '* 

= E + E = E + i - e 1(4 = e 10 (5 4) 

(e-^-Zp) 2 e-^-Zp (e'^-Zp 1 ) 2 e'^-Zp 1 l^-Zp 1 J (e'^-Zp 1 ) 2 

Then, combining equations <5.3> and J5.4I leads to a more symmetric representation of the linear approximation of the caustics; 

1 1 



2 _ e 3i * e 3i,t 
q (c = (4* - z P ) 2 + (e'^-Zp 1 ) 2 



(l-Zpe-*) 2 (l-z-'e- 



- 1 



(5.5) 



Therefore, up to linear order, the shapes of the caustics due to the perturbation from the planet lying at z p and the one at z p ' are identical. 
Note that z p and zT both have the same argument but that their norms are related inversely to each other. It is also straightforward to show 
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that this symmetry further i mplies the magnification degeneracy between these two planet positions through equation J3.13I . Compare this to 
the degeneracy identified bv lGaudi & GouldHl997l) as their second class of discrete degeneracies of planetary microlensing - "whether the 
planet lies closer to or farther from the star than does the position of the image that it is perturbing." Note that this degeneracy also extends 
to systems of any number of multiple planets since the superposition principle applies to the linear approximation. 

If one compares the actual central caustics of a planetary microlens to its linear approximation in equation 15.51 . one finds that deviations 
between them take place for smaller q with |z p | < 1 than |z p | > 1, and therefore that the z p -z p 1 degeneracy appears not to be strict for q > 10 2 . 
This is because the magnitude of the actual perturbation term (\qf\ ~ q) for |z p | < 1 grows faster than the one for |z p | > 1, which grows 
\qf\ ~ q\z P \ 1 f° r kol ~ 1. To find a better approximation of the caustics for close-in (|z p | < 1) planetary systems, one may need to include the 
second order effects of the perturbation (see Appendix[C}. However, if < |z p | <k 1, it is possible to find an alternative perturbative approach 
for which the leading perturbation term behaves as ~ q\z p \ 2 , and therefore the linear approximation that can be used for larger q. Similarly, 
one can devise a different description of the system in which the leading perturbation term behaves as ~ g|z p r 2 , which can be applied for 
wide-separation (|z p | » 1) planetary systems. Furthermore, one can also show that these two types of perturbative description of system, in 
fact, lead to a certain pair of LCIP, which may be understood as a natural extension of the z p -z p 1 degeneracy. 



5.2 extreme binary lens 

It was |Pominild ri999b) who noticed that the approximation of binary lenses with extreme separatio ns by c e rtain c lasses of perturbed point 
mass lens systems reveals the underlying connection between various binary lens systems. |Aibrow et alj 120021) tried explicit but rather 
limited calculations based on the perturbative approach to demonstrate the presence of the magnification degeneracy between two types of 
binary lens systems. Here, this magnification degeneracy is studied in a more general way as an archetypal example of the LCIP pair. 

Let us first think of the situation in which the lens equation is given by equation <B6t with two component masses. By series-expanding 
the deflection term caused by one of the masses at the location of the other mass; 

r = qi - qi = - q[ + qi y (f ~ hi )k _ _ g; + qi + y i2(z-i zi ) k 

where 1 ZI and l Z2 are the positions of the two mass points. The radius of the convergence for the infinite series in equation 15 -6i is given by 
|z - 4, 1 < |/ z , - 1 Z1 \. After some substitutions of symbols, 

w = qf\z - 4, ) ; to = -L (t - 1 Z1 - T^y) = q; U2 (( - 4, ) - -0w 2 - l m T\ (5.7) 
yli \ k 2 ~ 4i / 9i 

which is basically the choice of new coordinate origins [l tl for the image space and l Zi — <?2(4i ~ 4,) 1 f° r me source space] and Einstein 
ring radius of the mass at l Zi as the new unit of angular measurements, it is easy to see that lens equation 15.61 is a specific example of lens 
equation 13. 1 i for a point mass lens under perturbations, 



w £-1 



y*-i*\ (5.8) 



where 

k/2 

V . = qiqx = ®(7 - J r ( * +2) (5 9) 

Here, l Wk = q~ >/2 (l Zk - l zi ), that is, the lens location in this new coordinate system. In particular, 

»-f«*-^ » = Q^- (5 - 10) 

If lyol ^ L then, the circle \w\ - 1 lies well within the region of the convergence of the infinite series, and also the perturbative approach for 
the caustics in the previous section is valid. If one considers the linear approximation alone, the perturbation series may be truncated after the 
second term. Note that, here, the second order effect of the leading term in the series [ly^l ~ \l„ 2 - / H[ |~ 4 ] is higher order than the linear effect 
from the second term [|yi | ~ \l w% - Z W] |~ 3 ]. In effect, the lens equation 15.81 describes the system as a point mass lens lying in the gravitational 
tidal field due to masses that are well outside of the region around its Einstein ring. 

Next, let us consider the multipole expansion of the deflection terms of the lens equation JB6t at some point z c ; 

qi qi qi v & _ l ^ 12 v &a - 



qi qi q\ y ik, - z*) qi y 

^ ~ Z 7-1 v - I ~ Z 7-7.2-1 if- 7 \k J - 5 Zj 



z - 4, z - 4 2 z ~ z = to ( z ~ z ^ k z ~ Zc "to (z ~ Zc ^ k 

+ qi qih, + qiki - (qi + ?2)z c v> qiO Zi - z c )* +1 + ?2(4 2 - Zc) 



(7-7..) k+2 ' ( ^ 



Z-Zc (Z-Zc) tl (Z-Zc) 

Here, the radius of the convergence for the infinite series is given by \z - z c | > rnax(|4, - zd, I4 2 _ Zc\)- If one chooses z c to be the centre of 
mass of the system, 
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z c = L + t a = , (5.12) 

<?i + qi qi+ qi q\ + qi 

then the dipole term vanishes. Subsequently, substituting symbols via 

w = ( qi +q 2 r ll2 (z-Z c ); UJ = {q l +q 1 y lll ^- Zc \ (5.13) 

which is the choice of the centre of the mass z c being the coordinate origin and Einstein ring radius corresponding to the total mass of the 
system being the unit of angular measurements, one finds equation J5.ll> to be consistent with being another example of a point mass lens 
under perturbations (eg. 13.11 . 

0} = w Y Q 2 m w- {k+2 \ (5.14) 

where 

A gift, -zj + q 2 (l 2 -z c ) k q^+qil, 

Qlt ~ (? 1+92 )<™+' = q i+ q 2 (5 - 15) 

qigk - IJ + qjM, - \J_ q^-qit + fai - - k - [Z&^C-ftH - k 

(qi + qif'Kq! + qi) M (<?, + qi) M ' ^ (?t + <? 2 )* ' W 1 ' 

are basically multipole moments. Here, again l vi = (qi + q 2 ) ~ I/2 (' ;i - - z c ) is the lens location in the new coordinate system, but the relation 
to the original representation differs from the previous case because of the difference of the new coordinate system. The coefficient for the 
leading term (the quadrupole moment) and the following higher order moments can also be written as 

r( / H'2 _/ wi) . U4-2 k ~ ^4 — : tt — ('„,,-'„,,) - Qi,—, rr— ] — (hv 2 ~ hn) ■ (5.1/) 



~ (qi+qif ' ~" - (qi+qiY ^ ~ (<?i + qtf 

Again, if |<2 4 | <k 1, one finds not only that the infinite series converges at |w| = 1 but also that the perturbation series is small enough 
for the linear approximation of the caustics given in equation 13. 1 2t to be reasonably valid. Similar to the previous case, one finds that the 
second order effect of the leading term in the series [\Q^ ~ |/„, 2 - / W1 | 4 ] is a higher order effect than the linear effect of the second term 
[ISsI ~ \lw2 ~ lw, I 3 ] an d so one may truncate the perturbation series after the second term when considering of the linear approximation. 

Finally, if one compares equations 15 .8t and 15. 14i . one immediately discovers that the two become a strict LCIP pair if = Q 2 k+i 
(k = 1,2,.. .). However, as noted above for each case, one may truncate each series after its second term when one is considering only linear 
effects. Hence, the conditions that y = Q\ and j\ = (2 8 , in fact, suffice to define an LCIP pair from equations I5.8i and <5.14> . Following 
the usual convention of parameter definition for the binary lens, that is, d = (q t + q 2 )~ ll2 (l Z2 — l zi ) as the projected separation between two 
component in units of the Einstein ring corresponding to the total mass of the system and q = q 2 /q\ > as the mass ratio between the two 
components, one can rewrite the leading coefficients of the series 

ro = aW ri = a7^ (5 - 18) 
24 = ^; e 8 = a^f^, (5.i9) 

where the subscripts ',' and ' p ' are used to distinguish the parameters associated with the tidal approximation (for \d,\ » 1) and the multipole 
expansion (for \d p \ «: 1). Then, the condition for two system to be an LCIP pair becomes 



(1 + q,)(d,dpf = ?t (1 + ip) ; (1 + ^^ = 11*. ,5.20) 
iv 1 ~~ ip 

This condition can be used to find the system (q,, d t ) that is degenerate in the first order to the system (q p , d p )\ 



ip . ,_ 1 1 + ip « on 

(1-?p) 2 d P (l - ip + ip) 1 ' 2 

or the system (q p , d p ) to the system (q,, d,)\ 

l + 2 ? ,-(l+4 ? ,)" 2 (l+4 ?t y/2-l 1/1+4 ? ,\ I/2 

ip = 5 = 1 5 ; < { p = t\ TM — • (5 - 22) 

2?, iq t d,\ i + it I 

Algebraically, there is a second solution for the inversion of the equation <5.21> . However, the second solution is related to equation 15.221 
by ^ «-» q~ and d p «-» —d p , which results in the identical system because each relationship on its own corresponds to a 180° rotation of the 
system (with respect to the centre of mass) and so the combination of the two becomes the identity transformation. Note that, if q <k 1, then 
it » q p + 2q 2 , q p ~ qt - 2q 2 , and d,d p x 1 + (3/2)q. That is, this LCIP pair is reduced to the z p -z p ' degeneracy of planetary perturbations. 
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Figure 1. Linear approximation for the caustics of quadr-octu-pole lens given by eq. 15-231 . The real axis (the horizontal direction) is chosen such that both 
coefficients are real (i.e., (j>\ = cj> 2 = 0). The dots represent the coordinate centre. In the top panel, the shape change of the tetracuspoid caustics with varying 
.v = cijc\ at fixed c\ is represented. The bottom panels, on the other hand, show hexacuspoid (or tricuspoid if c\ = 0) caustics with fixed values of eq. 15.261 
- it is not the same as an area since the curves cross themselves. 



5.2.1 central caustic 

Detailed discussion regarding the analysis of the shape of the (central) caustics of extreme binary lens is available in the literature (e.g., 
lDominiklll999bl : iBozzalhoOfJ) . Here, the basic properties of these caustics are reexamined as an example of the simplest LCIP caustics, 
essentially controlled by two parameters. 

Suppose that a point mass lens is perturbed by a null-convergent perturbation of the form of / = —C\Z - c 2 z 2 . From equation 14. H . the 
linear approximation of the resulting caustic is 

& = 1 ( Cl e 3i * + 3c ie -'* + 2c 2 e 4 '* + Ac 2 ^*) = ^e* [e 3 *^'' + 3e- i( ^' >] + |c 2 |e* [e 4l( *-« + 2^*-^] , (5.23) 

where -2<pi and -3^ 2 are the arguments of c t and c 2 respectively (i.e., Ci = |ei|e~ 2 * and c 2 = |c 2 |e~ 3 *). Up to the same order, this 
caustic is identical to that caused by its LCIP, / = C[Z~ 3 + c 2 z~ 4 - In general, e'^^ 2 ' is not necessarily real. However, if one restricts to the 
approximation for extreme binary lenses, one finds that <f>\ = <f> 2 . That is, c t = y and c 2 /c t = 71/70 for the tidal approximation, or Ci = Q4 
and c 2 /ci = Qg/Q^, and therefore, from equations <5.18> and 15 .191 . both (/>[ and </> 2 are the argument of d (i.e., rf /\d | = e* = e*). Then, 

f c e~ ¥c = 2|ci| [cos 3 (0- <p c ) -isin 3 (0-^ f )] + \c 2 \ {4cos 2 (0 - c )cos[2((^ - cp c )] - 1 - 4isin 2 (0 - <j> c ) sin[2(^ - tj> c )]] , (5.24) 

where c/> c = (pi = <p 2 . If c 2 = (and the real axis is chosen such that <f> c = 0), the expression is reduced to the one derived bv lKovnenil987bll . 
To find the cusp points feq. l3.14t . 

— = — hicie- 2 '^ - c i e 2 " 1 ') + Slc^e-^ - c 2 e 3 "")] = e 1 * {3|c,| sin[2(0 - 0,)] + 8|c 2 | sin[3(0 - 0,)]) = 0. (5.25) 
d(f> 2i L v 1 \ - /1 

Here, if one still assumes that <p c = cj>\ = 4> 2 , one immediately finds that a cusp forms with (f> - <p c = tin, where n is integer. For these cusps, 
!fv[£.e~*] = (-l)"2|ci| + 3|c 2 | and 3[£ c e~*] = so that they lie along the axis with the angle of <f> c and are separated by 4|ci|. Here, similar 
to 3[z], the imaginary part of z, %[z] : C -» R is the real part of z (i.e., %[z] = (z + z)/2, 3[z] = i(z - Z)/2, z = %[z] + i3[z]). In addition, 
one can show that there is either one more solution of equation 15.251 in < (p - <p c < n if 4|c 2 | < |cj | or two solutions in the same interval if 
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4\c 2 \ > |cj |, and therefore that the caustic has four (4|c 2 | < IcJ) 2 or six cusps (4|c 2 | > \c t | > 0) unless c\ = 0. If C\ = and c 2 t 0, it becomes 
tricuspoid - three cusps - since the 0-period for eq. 15.231 then becomes n rather than 2n. However, hexacuspoid caustics generated by 
equation J5.23I with 4|c 2 | > |ci| are doubly-wound self-intersecting curves (Fig.0, which cannot be the caustics of any real pure binary lens 
(note that certain ternary lenses may be approximated by an octupole dominated perturbation series and thus consistent with hexacuspoid 
caustics although <f>\ and 2 may not be the same any more). Hence, the tetracuspoid condition may be regarded as one of the limit for the 
approximation of an extreme binary lens by a perturbation series. The area enclosed by a tetracuspoid caustic can be found by 3 

5 = k §' iMc = k X ~ ic % d0 = _ (l nlC[{2 + 47tN j ' (5 - 26) 

where the negative sign is caused by the orientation of the current 0-parametrization of the caustic. With the restriction that |c 2 /ci| <K 1/4, 
one may regard c\ as the main parameter controlling the size of the caustics while ,s = c 2 /c t is an asymmetry shape parameter for the caustics. 
In Fig.0 the change of the caustics shape with varying value of s at fixed c t is represented. The dimension is intentionally omitted since it 
can be linearly scaled with c\. In particular, the distance between two cusps on the horizontal direction is given by 4|ci| and is independent 
of s. For completeness, examples of the caustics with s > 0.25 are also shown in the bottom of Fig.Qalthough they are not applicable for the 
approximation of the caustics of any binary lens. 

The fact that the shape of the central caustics of extreme binary lenses is controlled not by the leading coefficient of the perturbation 
series but by the ratio of its two first coefficients also implies that, when the wing of the lightcurve is not well observed, there may exist a 
continuous degeneracy of the binary lens model running along the (almost) constant shape parameter. While the behaviour of the highly- 
magnified part of lightcurve when the source is near the central caustic can constrain the structure of the caustic well enough to determine 
the shape parameter, the constraints on the relative size of the caustics with respect to the Einstein ring require the precise determination of 
the overall time-scale. Since the magnification near the centre is essentially scale-free (A ~ c -f-> eq. l3.13> when the size of the caustic 
is sufficiently small compared to the distance between the caustic and the source, one may rescale the whole lightcurve by changing the 
blend fractions. Therefore, without the detailed constraints from the wing, the overall scale factor is left to be unconstrained and the model 
exhibits a continuous degeneracy with strongly correlated time-scale, peak magnification, blending, and the size parameter of the caustics 
(which is the leading coefficient of the perturbat ion series). One can find an archetypal demonstration of this degeneracy in the modelling of 
the lightcurve of MACHO 99-BLG-47 lAlbrow et all2 002V In Fig.|2| one can not only find an example of a two-fold LCIP degeneracy but 
also the continuous degeneracy with the constant shape parameter. 

It is notable that, as q — > 0, the shape parameter is mainly controlled by d alone (i.e., s ~ d p or s ~ d~ l ). In fact, if q <K 1, the 
approximation of the caustic by equation I5.5i can be valid, and one can immediately find that the shape of the caustic is completely (up to 
linear order) determined by the planetary position z p while the planetary mass ratio q becomes an overall linear scale. Moreover, Taylor-series 
expansion of equation 15.51 either for \z v \ <k 1: 

U4>) = qz p + | [3z*e-* + f^e 3i * + 4z 3 e -^ + 2 z 3 e 4i * + 0(\z p f)] , (5.27) 
or for |z p | » 1: 



Ut>) = I + I 



e 3i0 e 4i0 e~ 2i< * 

— + 3— + 2— + 4^- + 0(]z P r 4 ) 



(5.28) 

approaches equation 15.231 with the shape parameter given by s = z p (for |z p | <k 1) or s = zZ l (for |z p | » 1) as expected. 



5.2.2 planetary caustic 

Note that the condition for the series expansion I5.8i to be valid for |w| ~ 1 is that |Z W2 - / H ,| = (1 + q) [/2 \d\ » 1. That is, if q » 1, 
then the tidal approximation can be used even if \d\ ~ 1. However, for this case, one can easily notice that |yo| = (1 + q~ l y l \d\~ 2 ~ 1 
(eq. 15. 1 8i so that the lensing behaviour of these systems cannot be properly described by the perturbative analysis although the subsequent 
higher order terms may be ignored (or viewed as perturbations on a lens with linear shear). If one truncates the series after the leading 
term, one finds that the system modelled by the lens equation J5.8I is a point mass lens su bject to a constant external shear given by -yo, 4 
which is generally referred to as a Chang-Refsdal lens after Chang & Refsdall Jl979tll984h . While the lensing behaviour of Chang-Refsdal 
lenses is for the most part analytically tractable, its study is beyond the scope of the current monograph. Here, it is simply noted that, for 
planetary microlensing, provided that the planet lies sufficiently far from the Einstein ring, there exist two types of caustics; one that can be 



2 Strictly speaking, at 4|c 2 | = |ci|, the caustics develop a butterfly catastrophe, and therefore, it has three simple cusps plus one butterfly catastrophe. 

3 Note that the area two-form can be written as dx Kdy = (dz A rfz)/(2i) = d(zdz)/(2i). Hence, from the fundamental theorem of multivariative calculus, the 
area bounded by (z) can be found to be (2i) _1 / d(zdz) = (2i)- 1 jzdz. 

4 The negative sign follows from the definition given in Appendixltfl which conforms to the usual convention found in weak lensing literature - the direction 
of the shear to be the direction of the eigenvector associated with the positive eigenvalue of y. Assuming t e R, '-yo' is negative real (eg. 15. 81 so that, 
according to this definition, the direction of the external shear due to the companion at the the position of the planetary caustics is along the imaginary axis, 
that is, the shear runs perpendicular to the axis joining two masses. 
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Figure 2. Fig. 2 of lAlbrow et"afl <2002h along with additional dotted lines showing the curves of constant shape parameter s with the value for their best-fit 
models. For the extreme close binary, it corresponds to \s\ = ICs/Gd = (1 ~ + ?p)~ kpl while it is given by \s\ = lyi/yol = (1 + It) W" 1 for the 
extreme wide binary. (The original version © 2002 the American Astronomical Society). 



approximated by equation 15 .5t - the central caustic, and the other, which is basically the caustic of Chang-Refsdal lens perturbed by the 
higher order tidal-effect terms - the planetary caustic. In general, the lensing effects due to planetary caustics are rather smaller than those 
associated with the ce ntral caustics so that they usually appear as small additional signals on the main lensing event iGould & Loeblfl992L 
iGaudi & Gouloll997l) . 



6 DISCUSSION 

6.1 What about the critical curves? 

So far, the discussion on the LCIP pair has been focused on the identical (up to linear approximation) caustics and the degenerate (total) 
magnification. A natural question to follow is whether they are indistinguishable (at least up to the linear order) in every sense. The answer 
is actually negative as one can easily discover the difference in the critical curves and the individual image positions. For the monomial 
perturbation of f(z) = az", one finds, from equation 13. 161 . the corresponding critical curves to be; 

zA(/>) = e ,(4 jl - |n|a|cos[(n + 1)0 + , (6.1) 

where (p a is the argument of a, that is, a = |a|e 1(i " . It is straightforward to show that equation <6. It is not invariant under the LCIP transformation 
of n «-» — (n + 2) and a <-> —a (equivalently, (j> a «-» n-(f> a ), but rather the amplitude of the cosine linear deviation term changes to e(n + 2)\a\/2. 
In addition, while the azimuthal positions of the images for the LCIP pair are the same since equation J3. 101 is identical for them, one finds 
that the radial positional deviations given by 

ir=j[f + £ ^)](eR) (6.2) 

cannot be the same for the two different perturbations. However, it should be noted that the leading order descriptions of the critical curves 
and the image radial positions are dominated not by the first order but by the zeroth order terms, and they are still identical up to their leading 
zeroth order. More importantly, despite this difference, the actual observational distinction between the LCIP pair would be rather difficult 
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Figure 3. Geometrical determination of the azimuthal positions and magnifications of images with a linear approximation of the caustics. For this example, 
the two cusps lying along the horizontal direction are positive so that the tangent to a caustic point on the bottom half points upwards and vice versa. For a 
given shape the caustics, if one can draw a tangent to the caustics through a given source position (represented by a dot), there exists an image towards the 
direction (represented by an arrow) parallel to that tangent from the image plane centre. Note that the image plane centre is not necessarily coincident with 
the source plane centre nor with the caustics centre. The magnification of the image is proportional to the inverse of the distance between the source and the 
tangent point. 



since the critical curve is physically unobservable and two dominant observable characteristics of the individual images - the magnification 
and the azimuthal position - are indeed identical up to the linear order of the perturbations. 

To find the magnification and the azimuthal position of the individual images, one must solve equation <3.10> . which is a transcendental 
equation for <f>, and thus it may not be possible to find algebraic solutions although it is straightforward to find its roots graphically. Inciden- 
tally, it is also possible to find them in a purely geometrical manner provided that one can draw its caustics (eq. l3.12t . From equation 13. 141 . 
one finds that (d( c /d(f>) || e 1(4 , that is, the parameter <f> also determines the tangent direction to the caustics. This, together with the fact that 
equation 13 . 1 1 i is real when tf> is the solution of equation 13 . 101 . also implies that, if one can draw a tangent to the caustic through the given 
source position, the angle of that tangent line <p is a root of equation 13. 101 . Hence, there exists an image with its azimuthal position with 
respect to the image plane centre being <p and magnification given by (2/)~' where / is the distance between the source and the tangent point. 
Examples of this procedure for the perturbation given by f(z) = az or its LCIP f(z) = —af 1 is found in Fig.0 For definiteness, a is assumed 
to be real (horizontal direction) in this example, and therefore, the two cusps lying along the horizontal direction are positive. Finally, note 
that not all images are actually solutions of equation J3.10L However, missing images, if any, are typically higher order solutions that do 
not significantly contribute to the total magnification provided that the linear approximation is valid and dominant and that the source lies 
sufficiently close to the centre (and therefore the caustics). 



6.2 What is so special about the point-mass lens? 

The discussion until now has been limited to the case for which the unperturbed lens system is described by the point-mass lens equation M.W . 
However, the perturbation approach outlined in Section|3|can be applied to any lens system whose solutions are known. In fact, the procedure 
in Section|3|can be generalized to the case for which the unperturbed lens system is given by any circularly symmetric lens system that can 
form an Einstein ring. Suppose that the system is described by the lens equation 13 . 1 i with the unperturbed part given by 
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4G M(<r) 



(6.3) 



where r = \z\ is the radial coordinate of image, if/(r) is a circularly symmetric lensing potential, and M(< r) is the mass enclosed within a 
given angular radius r. If the equation 

AG 



;M« T) : 



(6.4) 



possesses a positive root r = ro, one can find the linear approximation for the tangential caustics of the system by series-expanding the lens 
equation at zo = r e'* and following procedure similar to Section l3~T1 In particular, one finds 



Zo 

e^/W) + e-^/'froe-*) + 2 



/foe" 1 *) 



,3i<6 



(6.5) 
(6.6) 

(6.7) 
(6.8) 

(6.9) 



J = (2-j3 ) 

(M = _£ 
r 2 

and 

Zc(#) = e i* _ f I 
r 2(2-^)1 

in place of equations 13 .91 . J3. 101 . <3. 111 . J3. 121 . and 13. 16i . Here, /Jo = (din M/d In r)\ r=ro is the local power index for the mass M at the 
projected angular radius r = r . In other words, the shape of (the linear approximation of) the caustics is completely determined by the 
perturbation alone once they are scaled by the angular Einstein ring radius (note that eg. 16.41 actually defines r to be the angular Einstein 
ring radius for the given circularly symmetric system), and the unperturbed mass distribution only affects the image magnification as a 
common multiplicative factor, and therefore, the magnification ratios between images are also completely independent of the unperturbed 
mass distribution. Moreover, it is easy to establish the same LCIP relation f(z) = aro(z/ro)" and f(z) = -aro(z/ro)~ <n+2) , or more generally, 
that the pair of perturbation potentials with an inverse symmetry with respect to Einstein ring radius ro [i.e., S\j/'(r) and 6if/ p (r) = Stf/'ir^r/r 2 )] 
has the same expression for the linear approximation of the tangential caustics. 



6.3 Are there LCIP pairs for convergent perturbations? 



As has been noted, the LCIP relation originates from a certain intrinsic symmetry of harmonic functions, that is, their azimuthal structures are 
invariant under the radial distance inversion. This implies that the two-fold correspondence LCIP does not exist for convergent perturbations, 
for which the perturbation potential is no longer harmonic. However, examination of equations I4.5i and I4.8i reveals more basic properties 
of the LCIP, that is, the azimuthal Fourier coefficients of the perturbation potentials at the angular Einstein ring radius are the same, which is 
indeed the generalization of LCIP to include convergent perturbations. 
Let us think of the potential 



Stff(r, 0) = ^ [e„(r) cos nip + g„(r) sin nip] = - ^ 



z" . r" 
—E„(r) H E„(r) 



perturbing a circularly symmetric lens system. Here, e„ 
4G M\ 



e„ + ig„ . Then, the lens equation is given by 



■■z 1- 



where 

r 2n+l 

Mr) ■ 



AelC 2 r 2 



- 6a, 



rl +2 dr 



^ 2r dr\r"l z"-' 2r dr "' 

n L 

m 



Bilr) 



Z" 

m-2 
[0 

„2n-l 



'o 7 _ * 



(r"s„) 



(6.10) 



(6.11) 



(6.12) 



(6.13) 



Here, the calculation utilizes the fact that r 1 = zz and consequently that 3=r = z(2r'T l = r{2z)~ l ■ That is, for any real function F(r), one can 
relate the complex derivative and the ordinary derivatives via 3=F = F'd,r = zF'(2rY x = rF'(2zY l ■ To find the expression for the linear 
approximation of the caustics, one series-expands equation <6. Ill and its Jacobian determinant at zo = r e 1(4 as in Section |3~T1 After some 
tedious but none the less straightforward algebra, one can find 



J = 2(1- K sfl ) 



(6.14) 
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(6.15) 



in place of equations 13.1H and 13. 121 . where k s $ is the (local) convergence of the unperturbed circularly symmetric lens at r = r (note that 
Po = 2/fs,o at r = r ), and the subscript V indicates the value of the function at r = r [e.g., f n $ = f„(r ) and g„ t0 = g n (r )\. Then, 



/»,o = flr ) - 

Sn,0 = Siifo) 

so that 



drlr'J 



°w,0 

ro 



E n,0 



-4r4<*0 



°n,0 



8n,0 



or 

ftjfc] 
9[fc] : 



-z- 



(« cos w0 cos + sin ncf> sin 0) + (n sin w0 cos - cos n<A sin <A) 



• Qn 

(;? cos rc0 sin - sin /10 cos 0) H j- (n sin n0 sin + cos n<f> cos 0) 



o 



(6.16) 
(6.17) 

(6.18) 

(6.19) 
(6.20) 



Note that the expressions do not involve the radial derivatives at all. In other words, up to the linear order, the caustics are completely 
specified by e„ = e„{r Q ) + iftX^o) or the azimuthal behaviour of the potential 6if/(r ,(p) exactly at r = r , and independent of its radial 
behaviour even in its immediate neighbourhood. In terms of the deflections due to the perturbations (instead of the perturbation potentials) 
given as in equation J6. 121 . this translates into the linear approximation of the caustics being a function of f H $ + g n $ alone, and thus the two 
perturbations producing the same linear approximation of the caustics if /„,o + gnfi is the same. The null-convergent LCIP pair is a special case 
of a pair of two perturbations given respectively by (J„,g„) = (c„, 0) and (J„,g„) = (0, c„) where the c„'s are (complex) constants. In terms of 
the potentials, they are the pair of s' n = ajjjr^f and s p „ = a n (ro/r) n where the a„'s are again some common constants. 



6.3.1 ellipticity, external shear, and quadrupole moment 



In many applications of gravitational lensing, the lensing convergence may be elliptically symmetric, that is, K(r,(f>) = K a [a(r, 0)] is the 
function of r and (/> only through the combination a = rp~'(sin 2 </> + p 1 cos 2 (f>) [/2 , where p is a constant. It is straightforward to show that this 
yields concentric ellipses with eccentricity of ] 1 -p 2 | 1/2 as iso-convergence contours. The analysis of these types of lensing systems is of great 
interest, but also very much more complicated than the circularly symmetric system. However, if p = 1, the problem can be greatly simplified 
by applying Fourier analysis. Let us consider the radial Taylor expansion of /c(r, 0) with fixed 0, that is, n(r, (/>) ss K(r, (p) + (dK/dr)\ r=f (r - f). 
Suppose that the expansion was at f = 2 1/2 a(l + p~ 2 )~ l/2 . Then the convergence at r = ap(sin 2 <p + p 2 cos 2 <pY 112 is 



K a {d) as K(f,(f))\ 1 + 



dlnK 



dlnr 



1-p 2 

cos 2(4 

1 + p 2 



- 1 



(6.21) 



Here, by noticing that K a (a) a nd (d\n/(/dhir )\ r= f do not involve any azimuthal dependence, equation 16.2 li can be used to derive the 
quadrupole approximation (c ,f . . iKovne J 1 987^1 of the elliptically symmetric convergence 



<c(r, <f>) 



1+p 2 



! 1 1-p 2 din*- 



2 1+p 2 d In r 



cos 20 , 



provided that |1 -p 2 | «: 1. 

Suppose that the lensing convergence is given by 

K(r, <f) = K s (r) { 1 + f cos[2(0 - <*,)]} . 

Then, with the canonical boundary condition, the potential can be found, using Green's function, to be 



(6.22) 



(6.23) 



f(r,<f>) = 2 



In r df fK s (r) + J dr ln r ~ \\~2 J d? ^fr'W cos [ 2 (<* ~ 4>e)1 + r 2 £ y ^K s (r) cos[2(<4 - . (6.24) 

If one defines radial functions pa{r) = K s rdf, which is basically a scaled mass associated with the circularly symmetric part of the 
convergence, and pi(r) = J^K s fdf, then the potential (up to an addictive constant) may also be written in (assuming tfi e is constant) 



i/r(r,<4) = 2 l A + 



2 r d f^l-L\\ {\fr 2 p\{r) + r 2 fdf 

Jr r 2 T Jo J r 



/«'i</) 



cos[2(<4- cf> e )], 



(6.25) 
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where r is some fiducial radius, i//q is a constant that may be formally defined by i//q = /Uo(r) In r + £ n' ( r ) m fdr. In addition, one can also 
find the deflection function 

a = 2drf = 2^P-z + ^ fdrrVH^e- 2 * - Z f"d* 4^e 2 *. (6.26) 
r Z Jo Jv r 

One can note that the potential given by equation 16.251 can be considered as a circularly symmetric potential perturbed by a small perturbing 
potential of the form of equation <6.10> with n = 2 provided that both of the integrals are small. Assuming that the circularly symmetric part 
of the potential in equation 16.251 allows to form an Einstein ring at r = ro [i.e., 2yUo(/o) = r^], the linear approximation of the tangential 
caustics is given by 



r e' A 

\2 

I " ar I r 

fo 



r= r»iL)*m. t ,z-_r^m. (6 . 28) 

Jo r \r J r J ro r \ r I r 



Here, the caustic is a tetracuspoid so that there are four images for a source position within the caustic. Equation J6.27L from a comparison 
to equations J5.231 . 15.241 . and <5.26> . also implies that the cross-section of quadruple images is (3n/2)(f„ + g^) 2 - Furthermore, the lensing 
behaviour of the source near the lens centre is essentially described by two constants /„* and g£. In particular, when the source and the centre 
of the lens are perfectly aligned (£ = 0), there are four images associated with the caustics, and their azimuthal coordinates are given by 
<f> = nn/2 + <f> e where n is an integer. It is easy to verify that equation J6.14I is real and that 



(r + 8r)€* = 



\ + (-1)" ^° g ° ■* 



J = (-D" +1 4(l - K sfl )(f* + g*) (6.29) 



2(1 - K s>0 ) 

so that the resulting images are an alternating-parity equal-magnification two-fold symmetric crucifix-form quartet with an axis ratio of 
1 - |(/ * - #g*)/(l - k s0 )\. For more general source positions, once the caustic is drawn using equation 16.271 . the azimuthal coordinates 
and the individual magnifications of the images can be found by following the same procedure as in Section l6~T1 using equation 16. 14i . In 
addition, since the form of the perturbation here is specified, it is also possible to find the radial coordinates of the images by the equivalent 
of equation I6.2i 

e -i * , r -n fr" 4 * + >~o \(f n *-0n) cos[2(0 - d> e )] + Of* +a*)sm\_2(d> -6.)]} 
Sr = -± [C + r [e'<^*V " ^^]} = ~ 1 ° ln , \ (630) 

for the given radial coordinate of the images (p. Algebraically, <p can be determined from the condition that Sr is real, that is, 3 e~ 1(4 ] + (/"„* + 
gT) sin[2(0 - f )] = 0. In other words, for a given source position, <f> is determined by /„* + which is in accordance with the geometric 
construction of <p since it is the same coefficient for the caustics (eq. l6.27| . 

Note the linear approximation of the caustics (eq. l6.27| is identical (after a proper rotation of coordinate) to those caused either by the 
null-convergent compact quadrupole perturbation with the quadmpol e moment of (/"„* + ^ *) r o or by me null-convergent tidal perturbation 
by a constant external shear of -(/ * + gX) (c.f. lwitt. Mao. & Schechte3 fl995). In fact, the examination of equation 16.261 suggests that 
those locally (near Einstein ring region) null-convergent perturbations can be caused by the quadrupole moment of the mass distribution, 
either entirely enclosed within the angular Einstein ring radius (/ oc z~ 3 ), or lying completely outside of the angular Einstein ring radius 
(/ oc z). In general, the lensing behaviour for the source near the lens centre when the convergence is given by equation 16.231 with |£ | <k 1 is 
basically equivalent to that of the system with a circularly symmetric convergence k s under the perturbation of / = ro[/o (z/ro}~ 3 — ^g(z/ro)], 
where f~ = / *e 2 '*' and gl = ^e 2 *. This can be generalized to the elliptically symmetric system that is additionally affected by external 
shear and/or compact quadrupole mass moment. 5 For this case, / ' = / *e 2 * + Q A and g' a = ggS 2 "^' + jo, 6 and they are not necessarily 
parallel to each other any more. Here, note that Q 4 = |<24|e 2l *e and y = |-y l e2l * r are complex numbers. (Strictly speaking, they are complex 
number representations of certain symmetric 2x2 traceless tensors. The absolute value of the complex number corresponds to the positive 
eigenvalue of the tensor and its argument is the same as twice the angle between the corresponding eigenvector and the real direction.) Since 
equations J6. 141 and J6. 1 51 are still valid regardless of /„ and g„ not being parallel, one finds, in the place of equation <6.27L - note the factor 
of two difference in the definition compared to eq. 16. 121 - that 

Ufi = "^Y- \fo +# t |[3e- i ^ +) + e 3i( ^ ) ] = 2r e^ |/J +^| [cos 3 (0 - + ) - isin 3 (0 - + )] (6.31) 
and also that 



5 In much of literature, a quadrupole lens and a lens with shear are used interchangeably without distinction. Although this practice has a certain justification 
since both the effects due to the constant external shear and the compact quadrupole moment are associated with cos 2<f> terms in the azimuthal part of the 
potential, they are always explicitly distinguished in this paper. 

6 The sign for the external shear is chosen to follow the usual convention in strong lensing literatures fe.g.. lHuterer et all2004lt . This choice leads the lensing 
shear due to the external shear term to be '-yo>' as defined in Appendix IbI 
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^ = 2(1^) K + r ° ^ * - 6 '° J} = 2(1^ ' ( } 

in place of equation 16.301 . Here, /J + g\ = |/ 1 + |e 2,fli+ and /J - ^ f = [/" 1 - |e 21 *- . As before, fr 6 R for that is the azimuthal 
coordinate of the image for the given source position. For f = 0, one finds the azimuthal coordinates of images <f> = nn/2 + + and also that 
2(1 - K s0 )Sr = (-l) n r [f^ cos[2(0 + - </>_)]• This implies that the resulting image quartets are identical for the set of and g\ if [/„' + g Q ] \ 
and [f^ - g\ || cos [2(0 + - 0_)]| are constant. For a more general source position f # 0, while one finds that the azimuthal coordinates and 
the magnification ratios of images are completely specified by + g\ , the projection that controls the radial coordinates depends on each 
azimuthal coordinate <p of the image, so that they are in general not common for the four images. However, if the source is sufficiently close 
to the centre, the projections are more or less parallel to one another, and therefore, there still exists a certain near degeneracy of the images 
concerning the choice of \f^\, |^ f |, and their relative orientation in the underlying lens model. 
If the convergence is given by a power-law ellipsoid 

*[«(r,*)] = -(-) = -(-) (psm0 + cos^) = - (-) ^— j ^l-_cos20j , (6.33) 

it can be approximated by a quadrupole convergence in a form of equation J6.23I with 

P(ao\**( 2p 2 \ (2 -« /2 1-p 2 / p\ 

Here, < /3 < 2. In addition, 

so that one finds ro = 2 1/2 ao(l + p ') 1,2 by solving 2po('o) = r 2 ,. Furthermore, pi(r) = ^p.o(r) since £ is now constant, and 2/c Si o = p. Note 
that 2fio(r) = r 2 , ^r 13 and (dlnno/dhir) = (dlnfii/dlnr) = p. Finally, 

70 Jo r Q \r ) r ? 2j 209 + 2) l+p 2 4 2 +/T 

= r dr pxa ^ = „ ^ = _ 3, = lyj 
J, r \r/ r 2 J, 2(J3-2) l+p 2 4 

and thus 

.,2 



/o* + ^0* = i , o2 T7~S ; /o* -2o =~ 1 , „2 tTo , m - ^ 6 - 38 ) 



l+p 2 2+/3' 70 *° l+p 2 2(2+/3)' 

Hence, for a singular isothermal ellipsoid (SIE) lens (2k s0 = P = 1), with an ellipticity parameter e = (1 -p 2 )/(l + p 2 ) = 1 -p, one finds that 
fa + So = e /3 an d/ * -,<r * = — «/6, and thus, the four-imag e cross-section is given by (3n/2)(f * + g^) 2 = (7t/6)e 2 , and the axis ratio for a 
crucifix-form quartet is 1 - \(f * — ^)/(l - k s ,o)\ = 1 - («/3) (Keeton, Kochanek, & Seliak 1997). 
If this SIE lens is subject to perturbations due to the constant external shear, one finds that 

ft = ^e 2 * = 1 (,e 2 * + 4| 7o |^) - ^e 2 '^ ; ^ = « e * + Ne 2*, = 1 + %o| ^ _ (639) 

By comparing this to another SIE lens system that is affected by the compact quadrupole moment, fJ = (e/12)e 2 * + |<24|e 2l * e and #1 = 
(c/4)e 2l(4 ', one can conclude that the lensing behaviours of two systems are identical up to the linear approximation provided that the external 
shear and the compact quadrupole moment are related by |yo| = 3| Q4I and e&wa-tr) = -1 and the two ellipticity parameters by 



1 + !!Z2l e 2iWr-«M 



(6.40) 



In other words, the lens system that can be modelled by a SIE lens with an external shear can also be modelled by a different SIE lens with 
a quadrupole moment, the magnitude of which is a third of that of the external shear, provided that the system's departure from the circular 
symmetry is small. This can be further generalized to systems that are subject to an external shear and a compact quadrupole moment at the 
same time to yield a certain degeneracy of the lens model regarding the ellipticity, the shear, and the quadrupole moment (as well as their 
relative orientations). 

If one models an image quartet using a SIE lens with a constant external shear but no quadrupole moment, then/ ' +gl = |/ 1 + ,? f |e 2 "* + = 
(e/3)e 2i,fc + lyole 2 ' ' - and/ + - 8 \ = |/ f - ^ f |e 2 '*- = -(c/6)e 2i(4 ' - |y |e 2i ^. Hence, the four-image cross-section is found to be 

Y l/o' +slf =7l [j + + e|yo|cos[2(^ - e )]J . (6.41) 

It is easy to see that, for given e and |y |, the cross-section is largest [= 7t(e + 3|yol) 2 /6] if <py = <Pe while it is at minimum [= 7t(c - 3|yol) 2 /6] 
when <p y - <p e = n/2. In addition, one can find that 
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(5) - Irol 2 = l/o f + sl\ {3l/o f + 4[/ f - al\ cos[2(0+ - 0_)]} ■ (6.42) 

That is to say, for fixed \f£ + and |/ 1 -^J] cos[2(0 + - $_)], which produce identical crucifix-form image quartets, there exists a degeneracy 
running along a hyperbolic path in (e, |yo|) space. This may be seen as a prototype for the 'cancellation branches' of the 'U' shape degenerate 



1. 1 1997) 



path found in fig. 8 of Keeton et al. 1 1997). In reality, it is highly unlikely that the source is perfectly aligned with the lens centre so that the 



actual degeneracy becomes rather more complex. (Furthermore, as the ellipticity and/or the shear get larger, the whole perturbative approach 
will start to break down.) Nevertheless, the basic idea that there exists a certain set of /„ + and / f - g\ that yields similar (or nearly 
degenerate) image configurations (with the proper rotation of frames) is thought to be more or less valid, and this can lead to the observed 
degeneracy of the shear and ellipticity combination. 
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APPENDIX A: COMPLEX ANALYTIC FUNCTION 

Let / = u(x, y) + iv(x, y) be a function of z = x + iy in the region R c C containing zo, where u and v are real- valued functions of real variables 
x and y. If there exist continuous partial derivatives at z = Zo and they satisfy a differential constraint in the form of coupled differential 
equations given by 

u ,x = v ,y ' u ,y + v ..v = 0> (Al) 

one can show that / is complex differentiable, that is, there exists a unique complex derivative 
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, ,. f(z)-f(zo) 
f (zo) = hm 

2^zo Z - Zo 

at z = Zo- (Here and in the following appendices, the subscripted comma notation for partial derivatives is used, that is, u x = d x u = du/dx.) 
Furthermore, if the condition given by equation IA1I is met for all points in R, then the function / is complex differentiable arbitrarily many 
times (i.e., infinitely complex differentiable) at every point in R. In general, a complex-valued complex-variable function / is referred to as 
an analytic function (or a holomorphic function) in R iff is complex differentiable at every point in R. Note that / is analytic in R if and only 
if it satisfies the condition given by equation JA1> in R, which is usually referred to as Cauchy-Riemann condition. If / = u + iv is an analytic 
function, Cauchy-Riemann condition JAU implies that 

"» + Uyy = (Uj)j + (Uy)j = (Vj,),.V + ("V^J, = Vjy ~ V„ = 0, 

V.V.V + V ,yy = ( V ,.v),Y + (V,y),y = (-",,),, + (",.v),y = ~U„ + = 0. 

That is, both component functions u and v satisfy the (two-dimensional) Laplace equation. In general, a function F : R 2 — > R 2 is referred to 
as a harmonic function if it satisfies Laplace equation V 2 F = (9 2 + 9 2 )F = F xx + F TV = 0. In other words, both the real and imaginary parts 
of any analytic function are harmonic. 

If one considers a complex conjugation z = x - iy of a complex variable, the function / = u + iv may be regarded as a function of z and 
z, by the relations x = (z + z)/2 and y = i(z - z)/2. Then, formally, one may consider a partial derivative 

(U,x + i«,y) + i( v ,x + IVy) u ,x ~ Vy . U , + V , 

fi = u= + iv, = (u x x, + u iy y,) + i(v x x, + v y y£) = — '■ ■ — = — + i- ' 



2 2 2 

That is, f, = if and only if / is analytic. In other words, an analytic function may be casually understood as a complex function of z 
alone. (Roughly speaking, / is analytic if it can be expressed by complex constants and elementary differentiable real functions with their 
arguments replaced by the complex variable z.) 

If / is analytic in R containing zo, there exists a convergent Taylor series expression of / at zo 

f(z) = V(z - zo)"^^ = /(zo) + /'(zoXz - zo) + f -^(z - zo) 2 + ^^(z - zo) 3 + ■ ■ ■ , (A2) 
n\ 2 6 

«=o 

where / ( "'(z) = d"f/dz"- A Taylor series (absolutely) converges within its radius of convergence. That is, there exists a positive real number r 
such that the series in equation <A2> converges if |z-Zol < r. A complex function / is said to have a pole of order n at zo if lim z ^ ;o [(z-Zo)"'/(z)] 
is unbounded for non-negative integers m < n and (z - Zo) m /(z) is analytic in a region containing zo for positive integers m > n. If z 6 R is a 
pole of order n of a function / that is analytic in R - jzoK then / has a series expression of the form of 



= (z zo)"/(z) ZoT „ f 1 *r _ _ zo) „, = y 1^ ^_ _ 

(z-Zo) *—i m) dz'" z =a, (m + n)\ dz m+n 

m=0 y m=~n 



(z-zo)'" 



which converges if |z - zol < r and z t Zo where r is the radius of convergence. Note that the series expansion of a complex function of the 
form of / = 2m=-oo a m(z - Zo)"' is generally referred to as Laurent series. In other words, if the a m 's are coefficients of the Laurent series 
expansion of / at a pole of order n, then a m = for negative integers m < -n and a_„ # 0. If / is analytic in C except at isolated poles, the 
radius of convergence for a Taylor series at a point zo that is not a pole is the same as the distance to the nearest pole while that for a Laurent 
series at a pole is also given by the distance to the next nearest pole. 



APPENDIX B: LENS EQUATION WITH COMPLEX NUMBERS 

In general, the lens equation for a given lensing potential if/ is given by 

y=x-V x f(x), (Bl) 

where y is the angular position of the source and x is the angular position of the image. In addition, the lensing potential is related to the 
mass column density per unit cross sectional area of lens E via Poisson's equation 

V^=^X, (B2) 

c l 

where Z)r = D^D^D^ 1 is the reduced distance 7 . Note that one may also rewrite equation <B2i using the surface mass density per unit solid 
angle Y, A = as in V 2 !/' = (8nGD~}c~ 2 )I.A- Here, also used is that DrD ic \ = Dj\ 

By choosing some axis as a real axis, one may rewrite the lens equation JB II using complex numbers; 

( = z-a(z,z)- (B3) 
7 If D s = D L + D LS , then D R [ = D" 1 + Djj. 
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where z = X\ + ix 2 , £ = y\ + iy 2 and o- = if/ A + iif/ 2 . Here if/ j = if/ V| and so on. The complex deflection function a is related to the lensing 
potential if/ by 

tfii + iil/i a 

fi.z = ^,1*1,2 + ^.2 x 2.z = — - = 2' 

(iA.il + i«A.2i) — i(iA,i2 + 1^,22) 0A,n + fe) + i(iA,2i - <A,i2) 1^2, 47tGZJ Rv , 

= + a -2 x 2,z = 7, = 7. = Tj^W = 2 S - ( B4 ) 

2 2 2 c z 

Here, note that x\ = (z + z)/2 and x 2 = i(z - z)/2. (The above equations also imply that 2F= = Fj + iF 2 = (VF)i + i(VF)2 and 4F-= = V 2 F 
for any real valued function F.) Sometimes, a z and » = are also referred to as the convergence k and the shear y, respectively, of the lensing 
potential. Note the differential relation connecting the convergence and the shear, k= = a-, = y-. The convergence k is always a non-negative 
real number, which is related to the mass column density by k = (4nGD R c~ 2 )Y, - some authors also define a critical mass column density 
S c = c 2 (4nGD R y l , and then k = S/S c . On the other hand, the shear y here is a complex number (i.e., it has two independent components); 

(if/.u + ii/f,2i) + i((fr,i2 + 1^,22) OA.n - ^,22) + iOA,2i + ^.12) 
y = a, = a A x u , + a a x 2rz = = 

Irl 2 = = = 7 (^.n - <A-22) 2 + (A 2 i2- 

It is notable, however, that the corresponding shear in real 2-dimensional notation does not behave as a vector under the coordinate transfor- 
mation (in particular, rotation of coordinate axes). Rather it is a symmetric tensor of a rank of two with a null contraction - a real symmetric 
2x2 traceless matrix, that is, 

/ 3t[y] 0[y] \ = l(tn-^2 2if/ A2 \ / cos 20 sin20 \ 
7 \ D[y] -%[y] J 2\ 2if/ M if/, 22 - if/,n j \ sin20 -cos20 /' 

where 2<p is the argument of 7, i.e., 7 = |y|e 2l( *. The two eigenvalues of y are ±\y\ and its determinant is -I7I 2 = -77. Note that the unit 
eigenvector associated with the positive eigenvalue I7I is (cos <f>, sin <f>) T whose direction defines the direction of the shear. On the other hand, 
the second unit eigenvector (- sin0, cos<£) r that is associated with the negative eigenvalue —\y\ is perpendicular to it. It is also notable that 
the lensing shear field should satisfy the differential constraint that y x {= k-=) e R, that is, 27R 12 = 71 u -71 22 , where y R = %[y] = \y\ cos 20 
and 71 = 3[y] = \y\ sin 20. This constraint is the direct result of the fact that the lensing deflection can be described by a scalar potential - 
i.e., 7 ~ = k-= = 4V 2 k = 4V 2 (2V 2 if/) = 8V 4 and V 4 is a scalar operator. 

In the case of the null convergence (k = 0), equation JB2> becomes Laplace equation so that if/ becomes a harmonic function. Then, 
one can find a complex analytic function if/ c (z) = if/(xi,x 2 ) + \ip(x\,x 2 ), where a function ip satisfies the Cauchy-Riemann condition for if/ c (z)- 
Then, it is easy to show that 

. ifr 1 -iif/ 2 if 1 - iip 2 1 . . _ 

</f c = >Pc,z = + i<P,z = 2 + 1 — 5 = 2 _ ^ - ^2 + ^,1) = a. 

which also indicates that, for the case of the null convergence, a becomes also a complex analytic function since it is a total complex derivative 
of an analytic function. 

If the lensing potential is given by 

4Gm , , , , 
D„\c- 



il/ r + if/ r 2Gm r - n 2Gm - 4Gm 

= = ^ _ + _ = _ _ = ^ _ 

the corresponding convergence is found to be 

4nGm „, . _ 1 m „, , 

/f = 7<S 2 (x - /) = 5 (x - I), 

Aeif 2 S c D rel D R 

where S 2 (x) is the 2-dimensional Dirac delta function and / is the lens position in real 2-dimensional notation. In other words, £ = mD^ 2 6 2 (x- 
/), so that the potential describes lensing by a point mass in in empty space lying at / (note that Aei^R = Di). Moreover, since Poisson's 
equation is linear in mass (density), the lensing by a finite number of multiple point masses can be described by the sum of individual 
potentials. Then, the deflection function is found to be 

— 4G \— ' m k 



— _ *U yi m k 

^~ D rsl c 2 ^- z _J zk - 

With this deflection function, the lens equation JB3> is reduced to a particularly simple (dimensionless) form 

qk =M> (B6) 

provided that every angular measurement is rescaled by the angular Einstein radius (eg. 12.21 corresponding to some mass M. Here, m k is the 
mass of k-th component located at the complex (angular) position L . While it is customary that the lens equation <B6> is rescaled by the 
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Einstein ring radius corresponding either to the total mass of system Q^ k q k = 1) - usually when several component masses are comparable - 
or to the most massive component (q { = 1) - usually when the mass of the system is dominated by a single component -, it is equally valid 
to choose any mass M as a unit of mass measurement as long as it is connected to the unit of angular measurement 9 by D ra c 2 8 2 = 4GM. 

The local (differential) behaviour of the lens mapping can be studied from the Jacobian matrix of the lens equation, which is basically 
the linear transformation that approximates the lens mapping locally. For the complex lens equation IB3I . its Jacobian matrix is 

Here, £ z = 1 — a x = 1 — K £ R and £j = -a= = -y. Since this matrix is Hermitian - the corresponding Jacobian matrix of the lens 
equation IB II with real 2-dimensional notation is real symmetric -, it has two real eigenvalues, which are A ± = 1 - (k ± |y|). Note that 
A + + A. = 2(1 - k) and A + - A. = -2\y\ < 0. Thus, \\ - A 2 = A\j\(k - 1) and so \A + \ < \A-\ if k < 1 and vice versa. Suppose that the principal 
argument of y = \y\e 2 " f ' is 2(f> e [0, 27t). Then, one finds that E + = e 1 * and £_ = iE + are the unit eigenvectors of Jacobian matrix, whose 
associated eigenvalues are A ± , respectively. Note that these two eigenvectors of the Jacobian matrix are also the eigenvectors of the shear 
tensor feq. lB5> . This is a natural consequence of the fact that J = (1 - k)\ - y, where J is the real 2-dimensional Jacobian matrix and I is a 2x2 
identity matrix. Finally, Jacobian determinant is invariant under the transformation of the lens equation IB II to the complex equation IB3I . 
that is, 

d£ A d^ = (dxj + idy 2 ) A (dy, - idy 2 ) = -2idyi A dy 2 = -2'\J dx x A dx 2 = -2iJ | dz ^ dz j A ( ^"^ j = A dz. 
Hence, 

J = det j ; = a + a_ = (i - K - lrl)(i - k + irl) = (i - k) 2 - Irl 2 . 



APPENDIX C: SECOND ORDER PERTURBATION CAUSTIC 

With the lens mapping given by equation 13. It , the linear perturbation approach becomes inadequate or at least insufficient as the magnitude 
of the perturbation grows. Here, the discussion in Section l3~T1 is extended to find the appropriate expression of the caustic up to second order 
of the perturbation. 

When the lens mapping is described by equation J 3 - 1 1 . its Jacobian determinant is found as equation J3.4I or 

1 I f a) f (I) \ 

J = l ~zk + A — + — H 2 / ( '7 (1) . (CD 



If one Taylor-expands equation JC1> at zo = e ,(4 and keeps up to the second order terms of Sz = z - Zo and e, 

J = e(g, + I,) - e 2 g lgl + [2 + 6(g 2 - 2f ,)] e*& + [2 + e(g 2 - 2 gl )] j*Sz - ^Sz 2 - 4SzSz - 3e 2l <*<5z 2 , (C2) 
where g„(<p) = e (»+iW/M( e i*). Also used ^ 

d:J = —^+0(<E); d- z d z J = —^?+0(e); d 2 J = —^+0(e). 

z z z z ~ z z 

Next, one inserts Sz = (ex\ + e 2 A / 2)e 1 * into equation iC2\ : 

J = e(4xi + gi + gi) + f 2 [4X2 - gili + (82 + h)Xi - 2 (8i + 8i)Xi ~ 1Q* 2 ] + 0(e 3 ). 

Here,^i wdxi are assumed to be real-valued. To find the expression for the second order approximation of the critical curve, one solves for 
J = + <9(e 3 ), that is, 

, ,s g\ + ii , ... 2 x\ - (si + h)xi + gih (gi + gi) 2 (gi + gi)(g 2 + h) giii 
xm = — — ; X2(<P) = - 4 = ^2~ + 16 + — • 

Then, the critical curve is found to be Zc(0) = Zq + Sz = e"*(l + exi + e 2 ^ 2 ), which, up to the linear approximation, recovers equation 13. 151 . 
The corresponding caustic can be found as the (second-order) image of Zc(<p) under the lens mapping 13. li . To do this, first, one similarly 
Taylor-expands equation 13. II at zo = e 1 * and keeps up to second order terms of e and Sz', 

Z = Sz+ e 2{4, Sz - e 3, ^f - ef - ef^Sz = e 1(4 (e _i *& + e^Sz - e^St - 6g„ - eg^Sz) , 

and replaces Sz by (exi + e 2 / ^ 2 )e 1 * and takes up to the quadratic terms of e; 

U4>) = ^ [e(2ATi - go) + ^(2X2 -x\~ fori)] , (C3) 

© 2004 RAS, MNRAS OOO.fTlETl 



SYMMETRY OF PERTURBATION LENSING POTENTIALS 2 1 




Figure 4. Various approximation of the central caustic of the system described by the lens equation; f = z - z~ 1 - q(z, - l p ) where g = 0. 1 and z p = 0.2. The 
horizontal direction is same as the real direction. The parameters are chosen such that the linear approximation of Section l5"T1 deviates noticeably from the true 
caustic although the perturbative approach does not fail completely. The curves shown in solid lines are (i) the linear order approximation in Section l5"Tl (ii) 
the second order approximation in AppendixlHl (iii) the linear caustic of the perturbation series of Section l5*2l keeping only the quadrupole term, and (iv) the 
li near caustic of the perturbation series of Section l5^2l keeping up to the octupole term. Also shown in dotted lines are the 'true' caustic found by the method 
o flWitlHl99d) . 

which is again reduced to equation <3.17t and consequently equation 13. 121 if one keeps only up to the linear term. One thing to note regarding 
the second order approximation, albeit obvious, is that the expression is no longer linear in the perturbation so that the superposition principle 
does not hold. That is, to properly account for all the cross terms, one always needs to work with all components of the perturbation together. 
An example of the second order perturbation caustic is shown in Fig.|4|along with alternative approximations for the same underlying caustic. 
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